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ABSTRACT 



O ■ 

^vq ■ The bulk flow, i.e. the dipole moment of the peculiar velocity field, is a sensitive 

probe of matter density fluctuations on very large scales. However, the peculiar ve- 
r- — i locity surveys for which the bulk flow has been calculated have non-uniform spatial 

■ distributions of tracers, so that the bulk flow estimated does not correspond to that of 

a simple volume such as a sphere. Thus bulk flow estimates are generally not strictly 

Q . comparable between surveys, even those whose effective depths are similar. In addi- 

S-) ' tion, the sparseness of typical surveys can lead to aliasing of small scale power into 

r/3 ' what is meant to be a probe of the largest scales. Here we introduce a new method 

of calculating bulk flow moments where velocities are weighted to give an optimal 
estimate of the bulk flow of an idealized survey, with the variance of the difference be- 
, tween the estimate and the actual flow being minimized. These "minimum variance" 

■ estimates can be designed to estimate the bulk flow on a particular scale with minimal 
t-H | sensitivity to small scale power, and are comparable between surveys. We compile all 

major peculiar velocity surveys and apply this new method to them. We find that 
most surveys we studied are highly consistent with each other. Taken together the 
data suggest that the bulk flow within a Gaussian window of radius 50 /i _1 Mpc is 
QV 407 ± 81 km s" 1 toward I = 287° ±9°, b = 8° ± 6°. The large-scale bulk motion is 

consistent with predictions from the local density field. This indicates that there are 
significant density fluctuations on very large scales. A flow of this amplitude on such 

■ a large scale is not expected in the WMAP5-normalized ACDM cosmology, for which 
| the predicted one-dimensional r.m.s. velocity is ~ 110 km s" 1 . The large amplitude of 

• i— i . the observed bulk flow favors the upper values of the WMAP5 fl m h 2 -as error-ellipse, 

■ but even the point at the top of the WMAP5 95% confidence ellipse predicts a bulk 

flow which is too low compared to that observed at > 98% confidence level. 

C3 \ 

Subject headings: cosmology: distance scales - cosmology: h is the Hubble constant in units of 10 km s~ Mpc~ ). In 
large scale structure of the universe - cosmology: observation the gravitational instabi l ity paradigm (iFeldman et alj|200ll : 
- cosmology: theory - galaxies: kinematics and dynamics - IScoccimarro et all l200ll : IVerde et alj|2002ft . this motion is 
galaxies: statistics due to the gravity of structures on larger scales. For a galaxy 

at po sition r, the peculiar velocity v is given by (|Peeblesl 
Il993t) 

1 INTRODUCTION ^0.55 , , ( r ' - r ) 

A long-standing question in cosmography is the origin of ^ 4ir J m ^ ' |r' — r| 3 ^ ^ 

the ~ 600 km/s peculiar velocity of the Local Group (LG) 

with respect to the Cosmic Microwave Background (CMB). where (5 m (r) = (p — ~p)/~p~, and p is the average density of 
The motion of the LG with respect to th e "Local Sheet" i n the Universe, f2 m is the matter density parameter, and we 
which it is embedded is only ~ 60 km/s (|Tullv et al.l l2008). have used 0.55 instead of 0.6 for the power of f2 m in the 
thus most of the LG's motion is due to structures on scales pre- factor to improve accuracy for models with dark energy 
larger than the Local Sheet, i.e. beyond 5 /i -1 Mpc (where (Linder 2005). 
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The issue of the LG's motion and that of other nearby 
galaxies has important cosmological and cosmographical im- 
plications. Specifically, as shown by Eq.[T] the peculiar veloc- 
ities of individual galaxies are sensitive to the matter power 
spectrum over a wide range of scales. Indeed, apart from th e 
Integrated Sachs- Wolfe (ISW) effect (|Sachs fc Wolfdll967h . 
peculiar velocities are the only probe of the matter density 
fluctuations on scales of ~ lOOfo -1 Mpc and clearly the only 
dynamical probe in the low-redshift universe. A given power 
spectrum predicts the r.m.s. of the components of a galaxy's 
peculiar motion. For models with more power, i.e. a higher 
normalization, one predicts a larger r.m.s. velocity. 

For a single galaxy, the contributions to its motion 
arise from a range of scales: from the local (~ 5 /i~ 1 Mpc) 
to the very large ( > lOO/i -1 Mpc) scales. One may reduce 
the effects of small scale density fluctuations by studying 
the peculiar velocity of a larger volume using a sample of 
peculiar velocity tracers such as galaxies, clust ers, or Type 
la sup ernovae. Beginning with the work of iRubin et al.l 
(1976), a number of such surveys h ave been undertaken 
over the last couple of decades (| Dressier et al.l 



Lvnden-Bell et al. 1988; Aaronson et al. 1989; Willick 



Court cau 1992; Ha n fc Moul d 1992; Mathc wson et al 



Colless et al.l 



2001 



20071 ; ISpringob et al.ll2007F ) 



Hudson et al 



1987 



199C 



19921; 



Lauer fc Postman! Il994l; iHudson et all Il997l; IWillick et all 



19971; iGiovanelli et al.l Il998al; 


Hudson et al.l 


1999]; 


Dale et al.l Il999albl; IWillickl 1 19991; 


Courteau et al.l 


2000|; 



The simplest statistic that can be derived from a sam- 
ple of peculiar velocities is the dipole moment of the sample, 
also known as its bulk flow. It was quickly realized that the 
bulk flow was closely related to the amplitude of fluctuations 
on l arge scales, and could be used to test cosmological m od- 
els (|Clutton-Brock fc PeebI3ll98ll ; IVittorio et al.lll986T ). At 
face value, however, the surveys cited above yield appar- 
ently conflicting results: the measured bulk flow ranges from 
to ~ 1000 km/s. Note, however, that many of the above- 
mentioned surveys are sparsely-sampled, and that while au- 
thors quote the bulk flow of the sample, this sample bulk 
flow is often mis-interpreted as the coherent bulk flow of the 
whole volume occupied by the survey. 

The issue of sparse sampling, small-scale aliasing and 
their effects on stat i stics such as the bulk flow were first 
analyz ed by iKaiserl (|l988T). and later 



dl995h and others Jjuszkiewicz et al 
Feldman et"al]|2003l:lSarkar et al.ll2007 



Watk ins fc Feldmanl 
20001; iHudsonl 120031; 



Watkins fc Feldmanl 



20071 ; iFeldman fc Watkins] 120081 ). These studies addressed 
the issue of comparing sparse surveys both to each other (to 
check for consistency between different sparse surveys) and 
the comparison of sparse peculiar velocity samples with ex- 
pectations from cosmological models. One lesson from this 
work is that both sparse sampling and aliasing present an 
important effect that must be accounted for in interpreting 
the results, particularly those from sparse surveys such as 
clusters or SNe. 

Bulk flow estimates are essentially weighted averages 
of the individual velocities in a survey. Previous work has 
focused on a weighting scheme that produces a maximum 
likelihood estimate (MLE) of the bulk flow of a survey, an 



estimate that minimizes the uncertainties due to measure- 
ment noise but does not make any correction for the survey 
geometry . Thus the MLE bulk flow is obviously depen- 
dent on a given survey's particular geometry and statistical 
properties. In this paper, we instead address the question 
of how peculiar velocity data can be used to estimate a 
more general statistic: the bulk flow of an ideal, densely- 
sampled survey with a given depth. Our approach will be to 
calculate optimal weights which produce the best possible 
estimate of this statistic. An approach related to this ques- 
tion is th at of Zaroubi and co llaborators, w ho used Wiener 
filtering (|Zaroubi et al.lll999T l and variants (|Zaroubil I2OO2T) 
to reconstruct the matter density field directly from pecu- 
liar velocities. That work, however, was focussed more on 
the mapping of the density field and the measurement of 
f3 = /(fi m )/6, wh ere b is the bias para meter, than on the 
bulk flow (but see lHoffman et aTl feoOl). In this paper, our 
aim is somewhat different: to construct dipole moments that 
probe the largest scales. 

The goal of this paper is to make the cleanest mea- 
surement of the large-scale bulk flow using the best peculiar 
velocity data available. We discuss the peculiar velocity sur- 
veys used in this analysis in Section [2] In Section [3] we de- 
scribe the construction of the velocity moments, the power 
spectrum model, and the optimal weighting scheme used to 
estimate bulk flow components free of small scale noise. In 
Section [3] we apply these optimal weights to the data. In 
Section O we assess whether the optimally- weighted bulk 
flow results from different surveys are mutually consistent. 
In Section [S] we discuss the cosmographic implications of 
our results. In Section [7] we compare the measured bulk flow 
with expectations from cosmological models. We discuss our 
resuts in Section [5] and conclude in Section 



2 DATA 

Here we analyze all of the recent peculiar velocity surveys. 
The datasets occasionally have outliers, and so it is neces- 
sary to remove them. Simply removing outliers with large 
CMB velocities might bias the resulting flow. Instead, we 
use predictions from the IRAS-PSCz density field to identify 
outliers, according to the following procedure. First, we com- 
pare the observe d peculiar velocity wi th the predicted pecu- 
liar velocity from lHudson et ail (|2004l ). adopting the param- 
eters of the B05 (/3 = 0.5) flow model used by iNeill et all 
(2007). This model allows for a small external bulk flow aris- 
ing from large scales, but provides a better predicted pecu- 
liar velocity within the PSCz volume, i.e. within a distance 
of 200 /i _1 Mpc. Peculiar velocities that deviate by more than 
3.5 a are rejected, where the uncertainty includes the dis- 
tance error and a thermal component of 150 km s _1 . For 
each sample, we also quote a characteristic depth defined as 
the mean weighted distance, where the weight is the inverse 
square of the peculiar velocity error. 

The final samples are as follows, listed in order of char- 
acteristic depth from nearest to most distant. 

• SBF: the su rface brightness fluctuation survey of 
iTonrv et al"1i|200ll ). We use the distances from their Table 1, 
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except where SBF galaxies are identified as group members, 
in which case we use the group peculiar velocities (their Ta- 
ble 4). After rejection of the outliers N4616, N4709 and ESO 
323-034, there are 69 field and 23 groups, with a character- 
istic depth of 17 ft - Mpc. 

• ENEAR: a survey of Fundame ntal Plane (FP) dis- 
tances to nearby earlv-tvpe galaxies (|da Costa et al 
iBernardi et~ai1l2002l ; IWegner et al.ll2003l ) . After the exclu- 
sion of 4 outliers, there are distances to 698 field galaxies or 
groups (Bernardi, priv. comm). For single galaxies, the typ- 
ical distance error is ~ 20%. The characteristic depth of the 
sample is 29 /i -1 Mpc. Note that unlike other samples con- 
sidered here, the se data are not corrected for inhomogeneous 
Malmquist bias (iHudsonll 1994a! ) . 

• SN ar e 103 Type la super novae distances from the com- 
pilation of iTonrv et al.l (|2003T ). limited to a distance of 150 
/i _1 Mpc. SN distances are typically precise to 8%. The char- 
acteristic dep th of the survey is 32 h~ Mpc. 

• SFI++ i|Springob et all l2007h . based on the Tully- 
Fisher (TF) relation, is the largest and densest peculiar ve- 
locity survey considered here. After rejection of 38 (1.4%) 
field and 10 (1.3%) group outliers, our sample consist of 
2675 field galaxies and 726 groups. For some analyses, we 
split this large sample into a field (SFI+- |-_f) and group 
(SFI++g) subsamples. The characteristic depth of SFI++ 
is 34 ft _1 Mpc. 

• SC l|Giovanelli et al.lll998bl ; iDale et al.|[l999ah is a TF- 
based survey of spiral galaxies in 70 clusters within 200 
/i -1 Mpc. The characteristic depth of the combined sample 
is 57 /i _1 Mpc. 

• SMAC (|Hudson et al.lll999l . 120041 ) is an all-sky Funda- 
mental Plane (FP) survey of 56 clusters. The characteristic 
depth of the survey is 65 fe~ 1 Mp c. 

• LP l|Lauer fc Postman|[l994l ; |Postman fc Lauerlll995l ) is 
a survey based on using brightest cluster galaxies as distance 
indicators. The survey consists of BCGs in 119 Abell clusters 
across the whole sky within a distance of 150 ft _1 Mpc. Here 
we obtain peculiar velocit ies using the metho d ology, but not 
the X-ray correction, of iHudson fe Ebelind (|l997l ). which 
makes a small correction to the error estimates of LP. The 
typical error per measurement is 19% and the characteristic 
depth of the survey is 84 fe~ 1 M pc. 

• EFAR (|Colless et al.l l200ll ) is a survey of 85 clusters 
and groups, based on the FP distance indicator. The EFAR 
survey was not intended to measure the dipole moment, but 
rather to examine peculiar velocities in two superclusters: 
Hercules-Corona Borealis and Perseus-Pisces-Cetus at a dis- 
tance of ~ 120 ft Mpc. As a result of this strategy, the 
coverage is far from all-sky. The characteristic depth is 93 
/i^ Mpc. 

• IWillickl l| 19991) is a Tully-Fisher based survey of 15 clus- 
ters with a characteristic depth of 111 ft" 1 Mpc. 

In addition to treating each of the above surveys in- 
dependently, it is also interesting to combine them into su- 
persets. The distance range spanned by the surveys is rather 
heterogeneous, however. Essentially the surveys fall into two 
categories: dense, relatively shallow surveys of nearby field 
galaxies or small groups (SBF, ENEAR, SFI++) and sparser 



but deeper surveys of clusters (EFAR, SC, SMAC, Willick). 
The SN sample straddles a range of depths, but is rather 
sparse, so we associate it with the latter category. However, 
the large numbers of objects in SFI++ dominate all samples, 
hence our superset labelled SHALLOW consists of SBF and 
ENEAR only, and we combine SFI field and group samples 
separately into a second shallow set labelled SFI++. The 
DEEP sample includes all other surveys, except for LP (see 
Section [S]). 

Finally, we also combine all surveys (except for LP) into 
a master catalogue labelled "COMPOSITE". The COM- 
POSITE catalogue has a characteristic depth of 33 /i _1 Mpc 
and is based on 4481 peculiar velocity measurements, mak- 
ing it the largest peculiar velocity catalogue studied to date. 



3 VELOCITY MOMENTS 

The statistics of individually measured galaxy or cluster pe- 
culiar velocities S n are not described well by linear theory 
due to the existence of nonlinear flows on small scales. This 
problem is typically solved by forming moments as linear 
combinations of peculiar velocities, u a — ^ w a , n S n , where 
W a ,n are a set of weights that specify the composition of the 
ath moment. For a proper choice of weights, and for a pe- 
culiar velocity survey that densely occupies a large volume 
of space, moments can be formed that are relatively insensi- 
tive to small scale motions and thus can be treated by linear 
theory; small scale motions are essentially averaged out in 
the summation. 

By far the most common moments used in the analysis 
of peculiar velocity surveys are the three components of the 
bulk flow vector. The bulk flow represents the net motion 
of the survey volume as a whole as traced by the galaxies 
occupying it. For an idealized survey, consisting of positions 
r„ and exact line-of-sight velocities v n for a set of N galax- 
ies or clusters, the bulk flow vector components Ui are just 
averages over the projections of the radial velocities onto the 
three coordinate axis directions, so that the weights for Ui 
are 

w it n = x, • r n /N. (2) 

For a more realistic survey, the measured line-of-sight 
velocity is assumed to have the form S n = v n + 5 n , where 
5 n is a drawn from a Gaussian distribution with variance 
+ a,. Here a n is the measurement error of the nth galaxy 
in the survey and cr* is the velocity noise, whic h accounts for 
smaller-scale flows not included in our model. iKaiserl (Il988l ) 
has shown that the weights for the maximum likelihood es- 
timate for the bulk flow components, which we will refer to 
as the MLE weights, are 

Wi.n = A/ ^ J. , "a ; ( 3 ) 

m 

where 

a \ " (x, ■r m )(x J -r m ) 

A »-l^ ^t^i • (4) 
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The statistics for velocity moments can be obtained di- 
rectly from the formulae for individual velocities obtained 
from linear theory. For example, the covariance matrix for a 
set of moments u a formed from the velocities S n of a survey 
is given by 



(UaUb) = ( 



^ W a ,nSn ] I ^ W b , m S m I } 
\ n / \ m / 



{S n S 



G n 



(5) 



The covariance matrix for the individual measured ve- 
locities Gnm = (SnSm) can be written in terms of the ve- 
locity field v(r) as 

Gnm = (v„V m ) + 5 nm (<j'l + <7„) 

= (f n • v(r n ) f m ■ v(r m )) + S nm (a 2 + a 2 n ). (6) 

In linear theory the first term can be expressed as an integral 
over the density power spectrum P{k), 

(f n ■ v(r„) f m ■ v(r m )> = 5=L J dk P{k)f mn {k), (7) 

where the function f m n(k) is the angle averaged window 
function 



f,„,,(k) = I — — (r„ • k) (f m ■ k) 



x exp (ikk ■ (r n - r m )) 



(8) 



Plugging Eq. ((6]) into Eq. (O and using equation 0, 
the covariance matrix of the moments reduces to two terms, 



(9) 



- Rab + R ab- 

The second term, called the "noise" term, is given by 

R ab = ^ w a>n w b , n (a n + cr») . (10) 

n 

The first term is given as an integral over the matter fluc- 
tuation power spectrum, P(k), 



(11) 



qi.i r°° 

Rab ~ 2^2 dk WUk)P(k), 

where the angle- averaged tensor window function is 

Waf>(fc) = ^ Wa,nW b , m J — — (f n • k f m ■ k) 

x exp (ik ■ (r„ - r m )) . 

(12) 

For the case a — b, Eq. (|12|l gives the angle-averaged 
window function for the moment u a . This window function 
tells us which scales contribute to the value of the moment. 
For moments that are measures of a component of the bulk 
flow, the window function should have the value of 1/3 at 
k = 0. This ensures that if the flow in the survey volume were 
uniform, i.e. all the power was at k ~ 0, the variance of the 



moment would correctly be 1/3 of the flow variance. Ideally, 
the window function for any useful moment should have a 
small amplitude at values of k corresponding to nonlinear 
scales; thus moments that measure the bulk flow components 
tend to have a peak at k = with the amplitude falling 
toward a plateau as k increases. 

For the MLE weights, the bulk flow moment window 
functions are determined by the spatial distribution of ob- 
jects as well as their associated velocity uncertainties. Given 
the fact that most surveys have relatively more objects at 
smaller distance, and that the measurement uncertainty in- 
creases rapidly with distance, the window functions found 
using the MLE weights tend to have broader peaks then one 
might naively expect given the depth of the survey, lead- 
ing to bulk flow moments that are sensitive to motions on 
somewhat smaller scales than the diameter of the survey. 

Since our goal here is to study motions on the largest 
scales possible, we require bulk flow moments whose win- 
dow functions have as narrow a peak as possible, also being 
small in amplitude outside the peak. Given that the mo- 
ments found using the MLE weights for a typical survey do 
not generally meet these criteria, we have formulated a new 
method for calculating weights for moments that essentially 
allow us to "design" the moment's window function, subject, 
of course, to the distribution and uncertainties of the data 
that is available. 

We begin by considering an idealized survey whose bulk 
flow components Ui have the desired window function. Here 
we will use an ideal survey which consists of a very large 
number of objects isotropically distributed with a Gaussian 
falloff in density, n(r) oc exp(— r 2 /2Rj), where Ri specifies 
the depth of the survey whose velocity is measured exactly. 
For this survey, Eq. <(2j) gives the weights for the bulk flow 
components, which will all have the same window function 
due to isotropy. Now, suppose that we have a galaxy or clus- 
ter survey consisting of positions r n and measured line-of- 
sight velocities Si with associated measurement errors a n . 
We would like to find the weights Wi, n that specify the three 
moments Ui = ^2 n Wi t „S n that minimize the average vari- 
ance, ((if, — Ui) 2 }. We will call these the minimum variance, 
or MV weights. The MV moments m calculated from these 
weights will then be the best estimate of the bulk flow of the 
ideal survey, if it were to exist, that can be obtained from 
the available data. We also expect that, within limits that 
will be described more fully below, the window functions of 
the m will match those of the ideal survey. 

In order to calculate the weights Wi, n , we first expand 
out the variance in terms of the weights, 

((Ui - Ui) 2 ) = Wi,nWi,m{S n S m ) + (Uj) 

-2^2w i>n (UiV n ), (13) 

n 

where we have used the fact that the measurement error 
included in S„ is uncorrelated with the bulk flow Ui. 

The next step would be to minimize this expression with 
respect to Wi. n ; however, as discussed above in order to be 
a proper measure of the bulk flow, the window function of 
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a moment must go to 1/3 as k — > 0. From Eq. (|12|l . we see 
that this requires the constraint 

lim Wii(k) = w itn Wi, m P nm = 1/3, (14) 

k — >0 ' * 



where 




We enforce this constraint using the Lagrange multiplier 
method; thus the quantity to be minimized with respect to 
Wi t „ becomes 



J2n,m W i.n w i,m{S n S m ) + (Uf) — 2 Wj, n (m Ui ) 




Taking the derivative with respect to Wi, n and setting 
it equal to zero gives 

^((5'„S m > + AP nm )w l , m = (S n Ui), (17) 

1)1 

which can be written in matrix form as 

(G + AP)w i = Q i , (18) 

where G nm = {S n Sm,) is the individual velocity covariance 
matrix, the components of Qi are Qi, n = (UiV n ), and w; 
is the TV-dimensional vector of weights specifying the ith 
moment. This is easily solved to give 

wi = (G + AP^Qi. (19) 

This equation, together with the constraint given in Eq. 
(|14[1 . allows us to solve for the weights in terms of the co- 
variance matrix G, the matrix P given in Eq. ((15}, and the 
n-dimensional vector Qi. Note that the MV weights can be 
positive or negative. If the ideal survey consists of TV' exact 
velocities v n i measured at positions r' n ,, then the elements 
of Qi can be written as 

JV' 

Qi,n = (UiV„) = i"!y(v»), (20) 

n' = l 

where w[ , = ■ r' n , /TV' as discussed above. The quan- 
tity (v n iv n ) can be calculated from Eq. (0 in terms of the 
positions r' n , and r n . In practice, we calculate Qi,„ by con- 
structing a simulated ideal survey where positions r' n , are 
selected at random to be isotropic and to have the density 
n(r) oc exp(— r 2 /2R 2 ). For the purposes of this study we 
have found that TV' = 10 4 points are sufficient for conver- 
gence of all relevant quantities. 

Note that the weights depend on the spectrum of mat- 
ter fluctuations, (see Eq. Q). Here we use t he power spec- 
trum model given bv lEisenstein fc Hul ljl998h . which explic- 
itly includes the effect of baryons. In this parametrization, 
P(k) oc k" s T 2 (k), where n a is the spectral index and the 
transfer function T(k) depends on the parameters h, Q m , 
Qb, the baryon density parameter, and <rg, the amplitude of 
matter density perturbations on the scale of 8 h~ 1 Mpc. 



4 RESULTS: THE MV WEIGHTS, WINDOW 
FUNCTIONS AND MOMENTS 

The minimum variance (MV) weights were calculated for 
the bulk flow component moments using the method de- 
scribed above for each of our catalogues. Here we will 
show results for two different ideal survey radii, a rela- 
tively shallow survey, Ri = 20 /i~ 1 Mpc, and a deep sur- 
vey, Ri — 50 /t~ 1 Mpc. For calcu l ating the weights, we as- 
sume the WMAP5 lDunklev et alj l|2008l ) central parameters 
n m = 0.258, Q b = 0.0441, a 8 = 0.796, h = 0.719, and the 
spectral index n s = 0.963, together with cr* = 150.0 km/s. 
We note that, for all but the sparsest surveys we consider, 
the values of the weights are insensitive to the specific power 
spectrum parameters used. 

One qualitative way to gauge how a moment con- 
structed in this way matches its ideal counterpart is to com- 
pare the window functions as calculated using Eq. (|12|l . Gen- 
erally, the larger and more geometrically similar a survey is 
to the ideal distribution, the better the match will be. For 
small surveys and/or those that have a very different spa- 
tial distribution than the ideal distribution the match can be 
quite poor. Measurement errors also play a large role in de- 
termining how good a match is obtained. Since the quantity 
((ui — Ui) 2 ) that is being minimized includes the noise term 
(uf), the optimal weighting scheme is a compromise between 
the need to have the moment's window function match the 
ideal window function and the need to keep the noise small 
by giving small weights to objects with large measurement 
errors. 

In Figures [T] - [4] we show the window functions of the 
MV bulk flow component moments. For comparison, we also 
include the ideal window functions as well as those for the 
MLE moments for each survey. As expected, the match be- 
tween the window functions for the MV moments and the 
ideal is best for the large surveys and those with small mea- 
surement error and similar distribution to the ideal sur- 
vey. For the sparse, noisy surveys, the window functions 
for the MV moments are not very different than those of 
the MLE moments, differing mostly in the amplitude of the 
tail of the window function for large k. There are some ex- 
ceptions: for example, for the Ri = 50fo -1 Mpc case, the 
SFI++ G and SHALLOW catalogue MV weights approxi- 
mate a R ~ 20h~ 1 Mpc Gaussian rather than the desired 
Ri = 50k- 1 Mpc Gaussian. The COMPOSITE catalogue 
window function is an excellent match on both 20 and 50 
Ti _1 Mpc scales. 

There are several ways to quantify how well the mo- 
ments constructed in this way should agree with their 
ideal counterparts. First, from Eq. (|13|) we can calculate 
{(m — Ui) 2 }. An alternative measure is to define the cor- 
relation coefficient (ui ■ Ui)/\ui\\Ui\. A value close to 1 indi- 
cates that the moment is an accurate measurement of the 
bulk flow of the idealized volume. While both of these mea- 
sures depend on an assumed power spectrum model, the cor- 
relation has the advantage of being dependent only on the 
shape of the power spectrum and not on its amplitude. The 
correlations between the MV moments and the ideal mo- 
ments given in Table [1] quantify how well they are expected 
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to agree, accounting for both agreement between their win- 
dow functions and measurement error. From this table it is 
clear that while our method works well for surveys which 
are large or have small errors, the moments it calculates for 
sparse, deep surveys that have large errors are not strongly 
correlated with the ideal moments that they are designed to 
measure. However, including these surveys in the composite 
catalogues does improve their correlations, particularly for 
the case of Ri = 50 h~ 1 Mpc. In particular, the DEEP cat- 
alogue does have a strong expected correlation (0.80) with 
the ideal survey with Ri = 50 /i~ 1 Mpc, even though each of 
its component surveys does not. As noted above, the corre- 
lation (0.91) of the COMPOSITE survey with the ideal sur- 
vey Ri = 50/i _1 Mpc is excellent, and the uncertainty due 
to the mismatch between ideal survey and actual weighted 
sampling is very small (~ 20 — 30 km s -1 ). 

In Table [1] we also show the magnitude of the MV 
moments, the expected deviation from the ideal moment, 
Ui — Ui) 2 }, given our power spectrum model. We also 
include the measurement noise in parenthesis to show how 
much of the deviation from the ideal moment is due to mea- 
surement noise and how much is due to differences in how 
the moments probe the window function. In Table [2] we 
show the Cartesian components of the MLE moments to- 
gether with the MV moments for both Ri — 20/i _1 Mpc and 
Ri — 50/i _1 Mpc. The Cartesian components as a function of 
the Gaussian window Ri for the composite surveys are also 
shown in Figure [5] The first thing to note from Tables [T] 
and [2] and Figure [5] is the remarkable consistency between 
the catalogues as to the value for the bulk flow components. 
This consistency will be explored in more detail below. 

It is important to note that for a given catalogue, the 
bulk flow moments calculated for different values of Ri do 
not have independent errors, and thus cannot be strictly 
compared. However, surprisingly, nearly all of the surveys we 
studied have a larger amplitude flow for Ri = 50 /i _1 Mpc 
than for 7?/ = 20 /i _1 Mpc (see Tables [Q and [2] and Fig- 
ure [5}. It is particularly compelling that this is the case for 
both the SFI++ and the DEEP catalogues, which are com- 
pletely independent and both of which have relatively small 
measurement errors on their bulk flow moments. A notable 
exception to this trend is the SC survey. This leads to the 
question of whether different surveys are mutually compat- 
ible, a topic addressed in the following section. 

5 CONSISTENCY OF SURVEY BULK FLOWS 
5.1 Method 

We now consider the consistency of the bulk flows measured 
by different surveys. While this has been done previously 
using the MLE moments, our MV moments have advan- 
tages. In particular, while MLE bulk flow moments depend 
on the details of a survey and are not necessarily comparable 
between different surveys, the MV moments have been de- 
signed to approximate the same window function regardless 
of the survey particulars. This means that the theoretical 
difference of the two moments will be as small as possible, 
making the comparison more rigorous. 



Following IWatkins fc Feldmanl (|l995l ), we quantify the 
agreement between two surveys, say survey A and B, by 
calculating the covariance matrix for the difference in value 
of the bulk flow moments uf and uf of the two surveys, 

D A-B n A B\i A B\\ 

Rij = {{Ui - Ui ){Uj - u 3 )) 

— R A J- R s _ R AB — R AB (01} 

ij ' ij ji ' \ J 

where the cross-terms are given by 

Ql.l /•oo 

R ^ = 2^J dk (W AB )%(k)P(k), (22) 
and 

(w AB )Uk) = J2 w t™f, m [ $ (*£•£ 

x exp (ik ■(r^-r„)). (23) 

Agreement of the moments from A and B can then be quan- 
tified by a x 2 f° r 3 degrees of freedom given by 

X 2 = - uf ){Rif B )-\uf - uf). (24) 

5.2 Results 

In order to quantify the agreement between the catalogues, 
we calculated the \ 2 f° r three degrees of freedom, as defined 
in Eq. 1241 for the difference in bulk flow of a given catalogue 
and that of the composite catalogue with the given cata- 
logue removed. Note also that since the covariance is model 
dependent, the results are given for the WMAP5 central 
parameters. The results are given in Table [3] 

From Table O we see that the MV- weighted bulk flow 
of LP disagrees with that of the COMPOSITE catalogue 
on both 20/i _1 Mpc and SO/i^Mpc scales. The level of dis- 
agreement on the latter scale, for erg = 0.796 corresponds to 
99% CL, although this would drop to 98% if erg were as high 
as 1. This is in agreemen t with previou s analyses based on 
MLE- weighted moments |Hudson|[2003). There are indepen - 
dent reasons to question the LP results: iHudson et al . (2004) 
compared, cluster by cluster, the distance to the brightest 
cluster galaxy derived by LP to that derived from the FP 
for other cluster members, and found that in a few cases, 
these distances differed significantly, in the sense that the 
LP BCG distance was too large (i.e. the BCG was fainter 
than expected). They found that all of the discrepant BCGs 
for which HST images were available showed strong evidence 
for dust. For these reasons, we have chosen to reject LP from 
the Composite catalogues. 

On the 50 ft -1 Mpc scale, the next most discrepant 
dataset is EFAR, which disagrees with the COMPOSITE 
catalogue at the 98% level, if erg = 0.796. The disagreement 
is more model-dependent, however, than is the case with LP. 
If, for example, erg were as high as 1, the disagreement would 
drop to only 92% C.L. Thus we choose to retain EFAR, as 
well as all of the other catalogues except LP, as part of 
the COMPOSITE catalogue. On the 20 h~ 1 Mpc scale, we 
note there is some tension between the SHALLOW cata- 
logue (dominated by ENEAR) and the SFlH — h catalogue. 
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Figure 1. The window functions of the bulk flow component for Rj = 20/i -1 Mpc for all of the catalogues we considered. The thick 
(thin) lines are the window functions for the MV (MLE) bulk flow components. The x,y,z-component are dash-dot, short dash, long dash 
lines respectively. The thick solid line is the ideal window function (since the ideal survey is isotropic, all component are the same). 



As noted above, on the largest scales, the DEEP catalogue 
and the SFI++ catalogue are in excellent agreement and 
both independently show a significant, large-scale flow. 



5.3 Comparison with kinetic Sunyaev-Zel'dovich 
bulk flow 

RecentljQ, iKashlinskv et all (|2008al fbh. have claimed a de- 
tection of the bulk flow from the dipole of the CMB observed 
behind clusters of galaxies, with amplitude (2.8 ± 0.7/iK) 



1 After this paper was submitted 
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Figure 2. The same as Fig. [T] for composite surveys 



and direction towards I = 283 ± 14° , b = 12 ± 14°. They in- 
terpret this as a dipole in the kinetic Sunyaev-Zel'dovich 
effect. The conversion from /iK to km/s has some system- 
atic uncertainty, but the authors interpret the bulk flow to 
be between 600 km s~ and 1000 km s _1 . 

Our bulk flow result is in excellent directional agreement 
(6°) with that found bv lKashlinskv et al.1 (|2008bh . The am- 
plitude of their flow (1000 km s -1 ) is considerably higher, 
but would be compatible if systematic and random errors 



reduced the Kashlinksy et al. result to ~ 400 — 500 km s . 
However, we note that their sample is very much deeper 
than ours. Whereas our signal arises from within a vol- 
ume of radius ~ 100/i -1 Mpc (z < 0.03), their signal is 
detected on much larger physical scales, with most of the 
signal arising from the shell in the range 0.04 < z < 0.2 
(120/1" 1 Mpc < r < 600/2T 1 Mpc) . 
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Figure 3. The same as Fig. □ for Ri = SOft-iMpc. 



6 IMPLICATIONS FOR COSMOGRAPHY 

Our robust measurement of a bulk flow of 407 ± 81 km/s 
toward I — 287° ± 9° , b — 8° ± 6° for a Gaussian window 
with Ri = 50/i _1 Mpc suggests that roughly 50% of the 
Local Group's motion is generated by structures beyond this 
depth. This is in good agreement with the value of 366 ± 125 
km s _1 toward I = 30 0°, b = 13° with i n a 6 /i~ 1 Mpc top- 
hat sphere, found by iHoffman et alj ([2001) based on the 
Mark III peculiar velocity catalogue. 



One way to locate the physical sources of peculiar ve- 
locities is through all-sky maps of densities of galaxies, with 
the assumption that one can use this map as a proxy for 
that of the dark matter mass, allowing for a bias factor, b, 
between galaxies and dark matter. The degenerate combi- 
nation of prefactors that scale Eq. [T]is then f3 = f(Q. ln )/b. 
Having mapped the density field, one can then predict the 
peculiar motion of the LG and other galaxies. By compar- 
ing these predictions with observed peculiar velocities, one 
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Figure 4. The same as Fig. |2] for Ri = SOft-iMpc. 



can solve for both /3 and for the bulk motion arising from 
sou rces at depths la rger than the galaxy survey. For exam- 
ple, iHudsonl Jl994bh mapped the density field of optically- 
selected galaxies within 80 /i _1 Mpc, a top-hat radius that 
is a close match to our 50 ft~ 1 Mpc Gaussian. By comparing 
the predicted peculiar velocities to the old Mark II pecu- 
liar velocity data set, he found (3 — 0.5 ± 0.06 and, more 
importantly, that the residual motion, arising from sources 
outside 80 /i _1 Mpc, was 405 ± 45 km s _1 toward I = 292°, 



b — 7°. This residual motion is remarkably consistent with 
the result found in this paper in scale, amplitude and di- 
rection (within 5deg), although it was derived from a com- 
pletely different data s et with differen t meth odology and 
assumptions. Similarly. iPike fc Hudson! (|2005l ). using a dif- 
ferent galaxy density field based on 2MASS photometry and 
published redshifts, found a slightly lower residual motion of 
271±104 km s^toward / = 300°, b = 15°. This result is also 
consistent with the Ri = 50/i _1 Mpc bulk flow found here. 
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Table 1. Bulk flow amplitudes |u| for the MLE-weighted case and for two choices of Ri for the MV weights. For the MV moments, the 
error in parenthesis is the noise error only. The quoted error includes both noise and the difference from the idealized survey geometry. Also 
shown are the correlation coefficients ( rjjrmr ) between ideal moments and their MLE estimates, where 1 indicates a perfect correlation. 
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Redshift surveys of IRAS-selected galaxies l|Strauss et al.l 
ll992l ; lRowan-Robinson et alj|200oh give very similar results 
for the residual flow beyond 60 ft -1 Mpc, provided Piras is 
~ 0.5, as found in d irect density-velocity comparisons (see 
iPike fc Hudsol]|2005l . for a summary table). These indepen- 
dent checks suggest that the large bulk flow motion is con- 
sistent with the absence of sufficiently massive attractors in 
the nearby (r < 80 /i~ 1 Mpc) Universe. 

Neither of the above mentioned galaxy surveys is deep 
enough to detect the physical source(s) of the 407 km s _1 
motion, which must lie beyond ~ 50/i _1 Mpc. The only all- 
sky galaxy survey t hat reaches much grea ter depths is the 
IRAS PSCz survey (|Branchini et al.lll999D . This shows lit- 
tle evidence of important attractors between 60 fa -1 Mpc and 
180 ft _1 Mpc, with the exception of the Shapley Concentra- 
tion, which is rela tively weak in IRAS. This issue was studied 
in more detail bv lHudson et al.l (|2004l ). who argued that, to 
explain the SMAC survey bulk motion, sources generating 
~ 200 km s _1 motion must be added to the IRAS-PSCz den- 
sity field. They attributed this to a combination of (i) sources 
in the Galactic Plane, (ii) sources beyond the 200 /i _1 Mpc 
depth of the PSCz, and (iii) to the undercounting of densest 
regions of the Shapley Concentration and the Horologium- 



Reticulum supercluster. The in terested reader is referred to 
Sec 5.4 of iHudson etafl |2004h for further details. 



On the other hand, if we assume that clusters of galax- 
ies trace the large-scale density field (perhaps with a higher 
biasin g factor) then one ca n use an all-sky survey of clus- 
ters l|Kocevski et alj 12004 ) as a predictor of the veloc- 
ity field. This survey suggests that, if clusters trace the 
mass, the Shapley Concentration and related very large- 
scale structures may play an important role, with as much 
as ~ 300 km s _1 arising from large scales, in approximate 
agreement with the large flow found here. 



Finally, it is worth noting that the large value of the 
residual motion implies that there are significant velocities 
generated by very-large scale structures and that this in 
turn has implications for the imp act of such structures as a 
"foreground" for ca l ibrating SNe JCoorav &: Caldwel|l2006l ; 
iHui fc Greenell2006l ; iNeill et al-lkoOTl ; lAbate fc Lahavll2008D 
and for its eff ects on the CMB such as the Integrated Sachs- 
Wolfe Effect |Fosalba fc Dore1l2007r ). 
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Table 2. Bulk flow vectors for the surveys (in Galactic Cartesian coordinates), for MLE- and MV weights . For the MV moments, the 
error in parenthesis is the noise error only. The quoted error includes both noise and the difference from the idealized survey geometry. 
Note that the MV Ri = 20h~ 1 Mpc and Rj = 50h~ 1 Mpc results for a given survey are not independent. 
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7 COMPARISON WITH COSMOLOGICAL 
MODELS 

As discussed above, the large-scale flow is directly sensitive 
to the large scales of the matter power spectrum, and so one 
can compare the observed value of the flow to that expected 
for a given cosmological model. Assuming that the Local 
Group does not occupy a special location in the Universe, 
one can calculate the expected mean flow and variance of 
bulk flow measurements taken at different locations. The 
former quantity is zero and the latter is the quantity of in- 
terest: it depends on the weights and sample geometry and 



is given by the covariance matrix Rij (Eq.0, which depends 
on the measurement noise and the power spectrum. This al- 
lows us to compare, in a frequentist sense, the observed bulk 
flow moments with the cosmological expectation. In partic- 
ular, we can calculate the \ 2 f° r three degrees of freedom 
corresponding to the 3 moments, as given by 

X 2 = ^UjR^Uj, (25) 

where i and j both go from 1 to 3 to specify the bulk flow 
components and the covariance matrix Rij is calculated as 
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Figure 5. The bulk flow of the composite catalogues as a function of Ri . Note that the data points are not independent, but rather are 
highly correlated. In the Galactic y direction (upper right panel), there is a consistent and robust flow exhibited by all catalogues that 
probe large scales, and this is reflected in the BF magnitude (lower right panel). 



described above for a specific set of values for the cosmo- 
logical p arameters. Here w e use t he ACDM power spectrum 
model of lEisenstein fc Hul l|l99ct ) and the WMAP5 central 
parameters as described above. While this statistic has been 
calculated previously for MLE bulk flow moments, the ad- 
vantage of the new MV moments is that, for the case of 
Ri — 50ft _1 Mpc, they have been designed to be sensitive 
only to scales of order 100/i -1 Mpc and larger. Thus we will 



be able to probe these scales without having to worry about 
the influence of smaller scales. Further, by isolating the very 
large scale motions, we will see that we will be able to put 
stronger constraints on power spectrum parameters. 

In Table|4]we show the expected r.m.s. bulk flow for the 
WMAP5 (tt m ,o- 8 ) parameters for the COMPOSITE cata- 
logue at a scale of Ri = 50ft" 1 Mpc. As can be seen by 
comparing Table [4] with the values in Table Q] the mea- 
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Table 3. \ 2 f° r 3 DoF for composite-surveys for Q m = 0.258. If the \ 2 value is greater than 7.8, the two surveys disagree at the 95% 
confidence level. 
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7.833 


10.796 


9.847 


DEEP-LP 


9.818 


8.351 


10.135 


8.964 


SFI++-LP 


8.025 


6.885 


11.541 


9.714 


COMPOSITE-LP 


8.405 


7.388 


12.344 


10.412 



sured and expected values differ significantly. Quantifying 
the disagreement, we find, for the WMAP5 central param- 
eters, find that \ 2 = 11-52. The probability of observing a 
bulk of flow this high an amplitude, in a Universe described 
by a WMAP5-normalized ACDM model, is only 0.9%. 



In order to assess the effect of the uncertainties in the 
WMAP5 parameters, we have explored further the multi- 
dimensional cosmological parameter space. Our covariance 
matrix Rij is dominated by the cosmic variance term (typ- 
ically ~ 100 km s _1 ) and not by the noise term (typically 
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Table 4. Expected bulk flows and observed y2 for 3 bulk flow moments (Rj = 50h~ 1 Wlpc), as a function of cosmological parameters 





ML 




BC 






Q m = 0.258 


cr 8 = 0.796 


n m = 0.262 


erg = 0.863 


Expected ID r.m.s. 


112 km/s 


121 km/s 


Survey 


X2 


P(> X2) 


X2 


P(> X2) 


SHALLOW 


1.54 


0.6731 


1.37 


0.7126 


DEEP 


7.54 


0.0565 


6.76 


0.0800 


SFI++ 


11.23 


0.0105 


9.92 


0.0193 


COMPOSITE 


11.52 


0.0092 


10.15 


0.0173 




WMAP5 + BA0 + SNA11 



99.737. (3 a) 



POSITE survey (Rj = 50k- 1 Mpc), given the observed flow, are 
shown by the black dashed lines at the 95%, 99% and 99.73% lev- 
els from top to bottom. Also shown ar e the W MAP 68% and 95% 
confidence limits from iDunklev et al.l ||2008|) (blue solid) well as 
for WMAP5+BAO+SN ijKomatsu et al.ll2008h (red dashed). The 
stars indicate the regions of the WMAP5 parameter space that 
maximize the bulk flow variance, and hence minimize the \ 2 - 



~ 40km s _1 ) which has a small effect when added in quadra- 
ture. Thus, to a good approximation, Rij should scale with 
the amplitude of the power spectrum, parametrized by erg. 
The cosmic variance term also depends on the /(fi m ) pref- 
actor and the power spectrum shape parameter, F, which 
on the large scales probed here is well approximated by 
r = flmh. Lower values of F lead to larger flows at fixed erg 
and /(fi m ). There is some cancellation of the f2 m -dependent 
terms, and so we have found that the flows depend on 
the combination fi m /i 2 . Fig. [6] shows the 2D parameter 
space of ftmh 2 and erg with h, Qb and n s are held fixed 
at their WMAP5 central values for this calculation. The 
large scale flow found here favours the corner with low Q m h 2 



and high a s . Also plotte d are the WMAP5 6 8% and 95% 
confidence regions from IDunklev et al.l (I2008TI. as well as 
the W MAP5+BAO+SN combination from iKomatsu et all 
(j2008T ). One can see that the entire WMAP5 95% CL pa- 
rameter space is excluded at better than 95% CL (2a). 
The WMAP5 "best case" (BC) par ameters, i.e. those that 
lie within the IDunklev et all l|2008h 95% CL regions, have 
0, m h 2 = 0.136 and as = 0.863 but the r.m.s. flow differs only 
slightly from the central WMAP5 value, and hence the \ 2 
value is very similar (see Table [4]). Essentially both WMAP5 
and flows are sensitive to the same large scales, and so they 
have the same parameter degeneracies. 

Another approach we took is maximum-likelihood anal- 
ysis. The cosmic variance in the bulk flow is a function of 
the cosmological parameters, and so one can ask which pa- 
rameters maximize the probability of generating a flow equal 
to that observed. Basically, given a set of power spectrum 
parameters 0, the likelihood is given by 



£(©) 



I.RI 1 / 2 



exp 



1 D-i 



(26) 



where i and j both go from 1 to 3 to specify the bulk flow 
components. Taking the natural logarithm of the above, and 
multiplying by —2, we see that this is very similar to the 
"frequentist" x 2 statistic (Eq. 125 p except for a term In \R\ 
that penalizes models with high cosmic variance. 

From Figure [6] we see that the bulk flow is more weakly 
dependent on Q, m h 2 than erg. This can be understood by 
considering the fact that increasing fi m increases the /(fi m ), 
but also increases F; these two changes act on the bulk flow 
in opposite directions and tend to cancel out. Given the 
stronger dependence on erg, we have chosen to plot the like- 
lihood of erg with other cosmological parameters fixed at the 
central WMAP5 values. The results are show in Fig. [7] The 
peak of the likelihood is at erg ~ 1.7, but the likelihood is 
very non-Gaussian: it has a sharp edge at low erg and has 
a tail extending to very large values. The 1-sided 95% and 
99% lower limits are 1.109 and 0.878, respectively. Of course 
these limits would be reduced if one placed priors on erg that 
excluded very large values. 
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95% lower limit 1.109 
99% lower limit 0.878 




Figure 7. Likelihood of the value of erg when the power spectrum 
shape is fixed at the WMAP5 central parameters. The dashed 
lines indicate the 99% and 95% lower limits. 



8 DISCUSSION 

There are several possible explanations for the discrepancy 
between observations and model that we have observed. 

First, it is possible that the large observed flow is the 
result of a systematic error in the data. However, the inde- 
pendence of the distance indicators (TF, FP and SN la) and 
methodology of the various surveys, as well as the agreement 
between different surveys makes this quite unlikely. Indeed, 
examination of Table[2]shows that no one survey is "pulling" 
the COMPOSITE bulk flow. Furthermore, we have shown 
that, on the Ri = 50h~ Mpc scale, the surveys are mutu- 
ally consistent within their random errors. Thus systematic 
errors affecting individual surveys must be small, and, since 
most systematic errors are expected to be independent, their 
net effect is smaller still. The only systematic that is likely 
to affect all of these surveys in the same way is a coher- 
ent (dipole-like) erro r in the foreground Galactic extinctions 
ijSchlegel et al.|[l998l ) used to correct magnitudes, and hence 
distances and peculiar velocities. This possibility, however, 
has be en tightly cons trained via extragalactic "colour" stan- 
dards l|Hudsonlll999l '). 

Second, there is the rather un-Copernican possibility 
that we happen to live in a rare local volume that has a 
statistically unlikely large bulk flow magnitude. The rareness 
would then be at the l-in-100 level. 

Finally, there is the possibility that the WMAP5- 
calibrated cosmological model underestimates the amplitude 
of large-scale fluctuations in the low- z Universe. In this con- 
text, it is interesting to compare our result to those of other 
independent low-redshift probes, which we consider from 
small to large scales. 



There are a number of small-scale measurements of erg 
(or f2 m °' 5 crg) from a variety of t echniq ues, some of which 
are summarized in iBond et alj l|2005l). The most r ecent 
weak-lensing result from the CFHTLS (|Fu et alj|200sh find 
tj 8 (fi m /0.25) - 64 = 0.785 ± 0.043 when all angular scales 
are analyzed, but cr 8 (n m /0.25) 53 = 0.837 ± 0.084 when 
only scales in the linear regime are studied. Some recent re- 
sults that sugg est higher va l ues of erg include the Lyman-a 
forest study of ISeliak et~aH ((2006), who find that that the 
power is enhanced over the expe ctation of WMAP3, and the 
study of iReichardt et alj (120081 ) which finds erg = OMtoM 
from secondary cluster SZ an isotropies in the CM B. Us- 
ing pairwise velocity statistics iFeldman et alj l|200St ) found 

„ — 1 1 q+0.22 
<7g - 1.13_ 23- 

It is possible to use peculiar velocities to probe matter 
density fluctuations on scales smaller than those probed by 
the bulk flow by comparing density-field predictions with 
observed peculiar velocities point-by-point. As discussed 
above, this yields an estimate of /3 = /(fi m )/b, which, when 
combined with an independe nt estimate of a s .oaU can be 
used to determine /(f2 m )erg. IPike fc Hudson! ([2005) com- 
bined their analysis of the 2MASS density field with previ- 
ous comparisons based on IRAS density fields. They found 
the degenerate combination erg (fi m /0.25) 55 = 0.88 ± 0.05. 
The quoted error is likely a slight underestimate, since not 
all studies are independent. Nevertheless, we note that the 
central value is somewhat higher than the WMAP5 value, 
although still lower th an our 95% lower limit. Although the 
IPike fc Hudsonl (|2005h measurement is based on peculiar ve- 
locities, the scale probed is quite different — in such studies, 
the large-scale bulk flow is essentially subtracted out, and 
so the scale is typically that of superclusters: ~ 20/t -1 Mpc. 
Similar statistical results based on peculiar velocities can 
be obtained from redshift-space disto rtions; for example, 
from SPSS LRGs, for fi m = 0.245, ICabre fc Gaztafiagal 
J2008h quote er g = .85 ± 0.06, in good agreement with 
IPike fc Hudsonl (|2005T I. but still in conflict with our formal 
95% lower limit on erg. In summary, from data on small to in- 
termediate scales, there are hints that erg may be higher than 
the WMAP5 best-fit value, but only by a modest amount 
which is not sufficient to comfortably explain the large-scale 
bulk flow. 

The bulk flow result result found here suggests signif- 
icant power on scales of fc < 0.02 h/Mpc. While the CMB 
probes such scales at very high redshift, there is only one 
other probe of the matter power spectrum at low redshift: 
the ISW effect. This p robe shows a signal that is stronger 
iGaztanaga et al J 12000 ), by a factor of 2.2 ± 0.6 l|Ho et all 
120081 ) than expected from the clustering of galaxies and 
quasars. It is interesting that our most-likely normalization 
erg is also a factor ~ 2 larger than the standard model. 

Fin ally, we consider the result of iKashlinskv et all 
(2008a. b). As noted by those authors, taken at face value, 
their flow amplitude and, more importantly, the scale over 
which this flow is observed (~ 6OO/1 -1 Mpc) is greatly in ex- 
cess of that expected in ACDM models. In this case, it seems 
impossible to generate cosmologically-consisent results sim- 
ply by tinkering with the parameters of ACDM; instead a 
wholesale revision of the model would be called for. 
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9 CONCLUSIONS 

We have calculated optimal "minimum variance" weights 
designed to measure bulk flows with minimal sensitivity to 
small scale power, and have applied these weights to a num- 
ber of recent peculiar velocity surveys. We find that all of 
the surveys we studied are c onsistent with each other , with 
the possible exception of the lLauer fc Postman! (|l994f ) BCG 
survey. Taken together the data suggest that the bulk flow 
within a Gaussian window of radius 50 /i _1 Mpc is 407 ± 81 
km s _1 toward I = 287° ±9°, b = 8° ± 6°. This motion is 
not due to nearby sources, such as the Great Attractor (at 
a distance of ~ 40/i _1 Mpc), but rather to sources at greater 
depths that have yet to be fully identified. 

A flow of this amplitude on such a large scale is not 
expected in the WMAP5-normalized ACDM cosmology. The 
observed bulk flow favors the upper values of the WMAP5 
S7 m /i 2 -(T8 error-ellipse, but even the point at the top of the 
WMAP5 95% confidence ellipse predicts a bulk flow which 
is too small compared to that observed at a confidence level 
> 98%. 

There are several possible explanations for the discrep- 
ancy we have observed. There is the possibility that we hap- 
pen to live in a volume with a statistically unlikely (< 2%) 
bulk flow magnitude. If this is the case, then the structures 
that cause this flow should be eventually identified as the 
depth and sky coverage of redshift surveys increase. Alter- 
natively, it is possible that the large observed flow is the 
result of a systematic error in the data, although the inde- 
pendence of the distance indicators (TF, FP and SN la) and 
methodology of the various surveys, as well as the agreement 
between different surveys makes this unlikely. 

The bulk flow in the nearby (d < 60k' 1 Mpc) Uni- 
verse is no longer noise-limited but rather cosmic variance 
limited, so that increasing the quantity of nearby pecu- 
liar velocity data will not alter the significance of this re- 
sult. At very large depths (d > lOO/i" 1 Mpc), however, 
the bulk flow measurement is still quite noisy. Future pe- 
culiar velocity surveys, such as the NOAO Fundamental 
Plane Survey (ISmith et all l2004h as well as nearby su- 
pernoyae surveys jlFilippenko et al.ll200ll ; Wqod-Vasev et al.l 
120041 ; iKeller et al.1 120071 ; iFrieman et al.ll2008h . are expected 
to yield more precise measurements of the amplitude of 
the bulk motion on these very large scales, and thus 
have the potential to strengthen the cosmological con- 
straints therefrom. In order to measure the bulk flow 
variance directly, one must measure the bulk flow in in- 
dependent (i.e. distant) regions. For the standard dis- 
tance estimators used in traditional peculiar velocity work, 
the errors grow proportional to distance and hence be- 
come infeasible at large distances. So other techniques, 
such as those based on th e kine t ic Sunyaev-Zel'dov i ch ef- 
fect |Sunvaev fc Zeldovichl Il972l; iRephaeli fe Lahavl Il99ll; 
Haehnelt fc Tegmarkl Il99rj; iRuhl et al.l |2004|; iKosowskv 



20061 ; IZhang et all 120081 ; iKashlinskv et al.1 l2008bl ) will be 



needed to access independent volumes. 

To reiterate, the results presented in this paper pose a 
challenge to the standard ACDM model with the WMAP5 
parameters. As this study shows, the implications to the 



standard scenario should be explored with as many inde- 
pendent cosmological observations as we can muster, with 
particular attention paid to clues from probes at low redshift 
and on the largest scales. 
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